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"Geometric scaling", i.e. the dependence of DIS cross-sections on the ratio Q/Qs, 
where Qg(Y) is the rapidity-dependent saturation scale, can be theoretically obtained 
from universal "traveling wave" solutions of the nonlinear Balitsky-Kovchegov (BK) 
QCD evolution equation at fixed coupling. We examine the similar mean-field predic- 
tions beyond leading-logarithmic order, including running QCD coupling. 



1 Motivation 



"Geometric scaling" (GS) is a striking empirical scaling property first observed in deep- 
inelastic (DIS) cross-sections . It consists in the dependence of j*p cross-sections on the 
ratio Q/QsiX), where logQs oc Y is the rapidity-dependent saturation scale. On a theo- 
retical ground, GS can be found as a consequence of saturation effects in QCD, when the 
density of gluons become large enough to impose unitarity constraints on the scattering 
amplitude. It has been shown [2] that the QCD evolution with a nonlinear term describing 
unitarity damping, the Balitsky-Kovchegov (BK) equation, leads to asymptotic "traveling 
wave" solutions exhibiting the GS property [2]. They are "universal" since they do not 
depend neither on the initial conditions nor on the precise form of the nonlinear damping 
terms. 

These results were mainly obtained at 
leading logarithmic order. In the present 
contribution, we describe how higher orders, 
in particular incorporating running QCD 
coupling, influence these predictions. Po- 
tential effects may be due to, e.g., next- 
to-leading (NLL) contributions to the evo- 
lution kernel, higher-order resummations 
schemes, observable dependence, infra-red 
rcgularization, position vs. momentum for- 
mulation. These aspects and the restaura- 
tion of universality at high enough rapidity 
Y has been discussed in Refs.[3], whose re- 
sults are here briefly described. The main 

difference with the fixed coupling prediction is a new kind of geometrical scaling with 
logQs oc VY, which appears to be as well verified by data [4] as the original GS prop- 
erty. 




Figure 1: Traveling waves 
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2 The Balitsky-Kovchegov equation with running coupling 

Before entering the discussion, let us introduce the traveling wave method in the case [2] 
where the running coupling has been introduced de facto in momentum space. One writes 

[bLog {k 2 /A 2 )] d Y T = {xll {-d Logk 2)}T - T 2 , (1) 

where T(k, Y) is the dipole-target amplitude in momentum space, xll the leading-log QCD 
kernel and [bLog (fc 2 /A 2 )] = a{k 2 ), the one-loop QCD running coupling. The asymptotic 
solutions of the BK equation can be obtained by recognizing the same structure [2] than the 
traveling wave equation u(t, x)—>u(t — vx) 

[x] d t u(t, x) = {d 2 + 1} u(t, x) — u 2 (t, x) , 

where the traveling- wavc/BK "dictionnary" is the following: 

Time = t ~ \/Y ; Space = x ~ logfc 2 ; Traveling wave u(t, x) = u(x — v c t) ~ T. 

Using the dictionnary, one thus recognizes the GS property u(x — v c t) ~ T(fc 2 /e 1 ' cV ^ 7 ), with 

a saturation scale Qs(Y) ~ e VcV ^ , where v c is the critical wave velocity determined [2] from 

the linear kernel xll- 

3 Traveling waves beyond leading QCD logs 

Let us introduce the general traveling-wave method for the extension beyond QCD leading 
logarithms. It consists in the following steps: 

• Solve the evolution equation restricted to linear terms in terms of a dispersion relation: 
u{t,x) ~ Jdrf e^-^*] 

• Find the critical (minimal) velocity v c = minw(7) = v{pj c ) which is selected by the 
nonlinear damping independently of its precise form. 

• Verify sharp enough initial conditions 70 > 7c, in order for the critical wave form to 
be selected. 

The mathematical properties of such obtained solutions ensure that the corresponding 
asymptotic solutions are "universal" that is independent from initial conditions, the nonlin- 
ear damping terms and from details of the linear kernel away from the critical values. Hence 
the traveling-wave method defines universality classes from which different equations admit 
the same asymptotic solutions. One caution is that the range of asymptotics may depend 
on the singularity structure of the kernel. This may have a phcnomcnological impact on 
the possibility of using these solutions in the available experimental range of rapidities. The 
saddle-point behaves as uj s ~ Y~ x / 2 except near singularities in 7 of the kernel. 

In order to illustrate the method, let us consider the general form of the NLL-extcndcd 
BK equation, replacing in equation |T|), xll — * Xnll (—dL,dy) , where L = logfc 2 . 
Introducing the function 

rl 

x (l, UJ )= &i Xnll{~i',w) , 

the linear solution reads 

-lL + -J- (2X( 7 , u„) - w s X( 7 , w,)) , 



T(L,Y) = J |£r ( 7 ) exp 
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where a dot means and uj s is given by the saddle-point equation 

Ybuj 2 s - X(j,UJ s ) +LU s X(j,LU s ) = . (2) 

From the solutions of the saddle-point equation, one can infer [3]: 

• For generic kernels beyond leading logs: The kernels may contain singularities up to 
triple poles due to the NLL contribution. By integration, new single and double-pole 
singularities appear in X at next leading order. The universality class is still the same 
but subasymptotics corrections may be large, and thus the critical wave solutions 
delayed to very large energies. 

• For Renormalization- Group improved kernels |5j?: The behaviour of the kernels xnll 
near the singularities are simple poles. This leads only to mild logarithmic singularities 
in the function X(-f, u>). The net result [3] is that one finds the same universality class as 
the equation {!]), since the oj dependence in X can be neglected and thus xnll^Xll- 

The same approach has been followed for recent QCD formulations of the Balitsky-Kovchegov 
equation with running coupling constant obtained from quark-loop calculation [6] . It leads 
to the same conclusion (with the same warning about eventual kernel singularities): the 
universality class for the BK equations with running coupling is the one defined by Eq. (fT| . 



4 Geometric Scaling in \/Y . 
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On a phcnomcnological ground, the main property of solutions corresponding to the univer- 
sality class of Eq.([T]) is the traveling-wave form u(x,t) ~ u(x — v c t) = u(fc 2 /e" cV ^ 7 ) in the 
asymptotic regime. Assuming a simple re- 
lation between that amplitude and the "f*p 
cross-section, one is led to look for geomet- 
ric scaling of the form cr 7 * p ~ (Q 2 /e Vc ^) 1 
with v c = est. In Fig|2j one displays the 
corresponding data plot |4] . The validity of 
the scaling property has been quantified us- 
ing the 'Quality Factor" QF method, which 
allows to determine the adequacy of a given 
scaling hypothesis with data independently 
of the form of the scaling curve [3] . 

One may also use the QF method to 
evaluate the scheme dependence of the sub- 
asymptotic, nonuniversal terms in the the- 
oretical formulae. In this case, the geo- 
metric scaling prediction is considered in a 

"strong" version, namely, with the critical parameters (such as v c ) fixed apriori by the the- 
ory. In Figj3] one displays the QF for geometric scaling for different NLL schemes. The 
top QF is larger than .1 which ensures a good GS property (similar than figE}. Depending 
on the resummation scheme (S3, S^.CCS, see [5]), |lo| gives the typical strength of the 
non-universal terms. The S4 scheme seems to reach GS sooner (at smaller |io|)- 
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Figure 2: Geometric Scaling in s/Y [4] 
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Let us give the main results of our analysis: 

• Mean-field saturation beyond leading- 
logs: The modified Balitsky-Kovchegov 
equations including running coupling and 
higher-order QCD corrections to the linear 
kernel asymptotically converge to the same 
traveling-wave solution. 

• Characterisation of the universality 
class: The universality class of these solu- 
tions is the BK equation with the leading 
logarithmic BFKL kernel supplemented by 
a factorizcd running coupling whose scale is 
given by the gluon transverse momentum. 
Higher order contributions to the kernel will 
affect the subasymptotic behaviour. 

• Higher- order effects in the kernel: The 

renormalization-group improved kernels are expected to improve the convergence towards 
the universal behaviour, spurious singularities being canceled. 

• Geometric Scaling: Geometric scaling in s/Y is a generic prediction of the universality 
class of the BK equation with running coupling, ft is well borne out by actual data, using 
the "Quality Factor" method [4] to quantify the validity of the scaling hypothesis without 
assuming the scaling curve a priori. 

• Nonuniversal terms: When using the theoretical "critical" parameters geometrical 
scaling is verified but requires the introduction of scheme-dependent subasymptotic 

Prospects of the present studies are interesting. On the theoretical side, it would be 
fruitful to investigate the universality properties of QCD equations beyond the mean-field 
approximation. On the phenomenological side, the problem is still not settled to know 
whether there is a slow drift towards the universal solutions or whether it exists subasymp- 
totic traveling wave structures, as mathematically [7] or numerically [8] motivated. 



Figure 3: NLL Quality Factors [4] 
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